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EXACT SOLUTIONS TO THE COMPRESSIBLE 
NAVIER-STOKES EQUATIONS WITH THE CORIOLIS AND 

FRICTION TERMS 

ANASTASYA KORSHUNOVA 



Abstract. We consider special solution to the 3D Navier-Stokes system with 
^^, ^ and without the Coriolis force and dry friction and find the respective initial 

^sj . data implying a finite time gradient catastrophe. The paper can be considered 

as extension of the results [I] . 

"a 

' O ■ We consider the following gas-dynamic like system: 

i 

^' (1) §^ + (v,V)p + /,(V,v)=0, 

3^; (2) p(^ + (v,V)v) = -Vp + p(Lv + AiAv) 



(3) ^ + (v,Vp)+7pdivv = 0, 

CO ■ where /5(i,x), p(i,x), v(i,x) — (wi,W2,V3) are density, pressure and velocity vector, 

. ■ respectively, fj, is the dynamic viscosity coefficient and 

o 

oo 

O ■ f-m -I o^ 

V 0; 

C^ ' 771 = const > is the friction coefficient, / ~ const is the Coriolis parameter, 

xgM", t> 0. 

Exact solutions of the system m)-© has been an area of intensive research 
activity in the last decades (e.g. [5], [3], [5], [5], [5]). Examples of exact solutions 
with a special initial distribution of the tangential component are given in [7] . 

Below we consider the solution to I!])-® in several particular cases. 
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1. Case I — 0, m — (without the Coriolis force and dry friction) 
Let us consider the velocity field and the density in the following form: 



(4) 



a{t)x 

my 

.W{x,t), 



p{t,x) = n{t) + ^{a{t)+P{t)>- 



Here a{t) + (5{t) is the divergency of velocity field. 
Firstly we are going to define the functions p{x,t), a{t) and P{t). 
For this class of solution the conversation of mass ([l]) and vorticity conversation 
equation 



(5) 



(9v 

Vx[p(— + (v,V)v)-/.Av]==0, 



which follows from ^, gives 



(6) a + ^(d + P)x + a^x{a + /3) + [a + ^(a + I3)x]{a + /?), 



[a+^(a + /?)a;] 



d'^W d'^W 



dxdt dx'^ 



^(a + P) 



dW 



ax- 



dW 
dx 



(7) 



dW 

dx 



a + TT(a + /3)a; 



dx^ 



(a + /3)(/3 + /32)=0. 
In [T] it was shown that a{t), [3{t) and ait) satisfy the following system: 



(8) 



m = 



t + B 



-a{t) ^ c{t) + m, 



(9) 



c{t) 



exp I — J a{T)dT j 

t 71 N 

C + / exp J a{T')dT' ) dr 
Vo 



where a{t) — exp J c{T)dT j and the constants B and C can be found from the 

initial data, namely: B = (/3(0))-\ C = (c(0))-i = -(a(0) + (3(0))-^ 

From ((8]) and ^ we find functions a{t), j3{t) and c{t). Substituting —a{t) for 

c{t) -\ — in ([52]) . we obtain: 

t -\- B 
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{t + B)exp (j c{T)dT 
(10) c(i) = - ^° 



C-J{t + B)exp ( / c{T')dT' ) dT 
Vo 



Let us denote p{t) = C — /(r + B)exp I / c{T')dT' j dr. It is easy to see that 

Vo 



c{t) = In 



p'it) y 



Thus, ([TOD impUes: 



(11) 



i + B 



p'it) Ci 



t + B p{t) ' 
therefore p^(i) = Ci(i^ + 2i3t + C2). Taking into account (3), we obtain: 

p'{t) C^{t + B) t + B 

c{t) - - 



p{t) p^{t) t^ + 2Bt + C2' 



C 
where Ci = - — , C2 = -BC. 

Therefore, we can find aU functions: 



(12) «(o = -r^ + T^^7^^, m ' 



t + B t^ + 2Bt + C2' t + B' 



(13) a(t) 



C. 



1/2 
2 



V<^ + 2Bt + C2 ' 



il/2 



/Oo„^_^^_ _ C2 , Po ^ + -8 



(14) «(a;. f) = cr(t) - ^c(i)a; , , 

^ ^ '^^ ' ^ ^ ' U ^ ' VWT2Bt + C^ U t^ + 2Bt + C2 

Having the analytic form of the solution, we can find initial data implying an 

unbounded increasing of its derivative in a finite time (gradient catastrophe) . It is 

clear that it is sufficient to find when the f{t) = — {t + B){t'^ + 2Bt — BC) has the 

positive zeroes. We have the following results: 

1) if B > 0, C < (/3(0) > 0, a(0) > -/3(0)) then f{t) 7^ 0, Vi > 0; 

2) B > 0, C > (/5(0) > 0, a(0) > -/3(0)) then /(t) has one positive zero: 



1 / / "(0) 



3) if B < 0, C < or C > -B (/3(0) < 0, a(0) > 0, a(0) ^ /3(0))then f{t) has 
one positive zero t = — _B; 
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4) if C = 0, S < (/3(0) < 0, a(0) = -/3(0)), then f{t) = (t + B){t^ + 2Bt). 
Thus, f{t) = in the points t = Ti = -B and t = Tj = -25; 

5) if B < 0, < C < -B (/3(0) < 0, a(0) < 0), then f{t) turns into zero in three 
points: 



i^T,^-5-v/^^T^.J^(^ 



a(0) _^^ 



(0)+/3(0) 



i = T3 = -5+v/^^T^=-^(^ 



«(0) , / 



a(0)+/3(0) 
It is obviously that Ti < T2 < T3. 

Now we are ready to find the third component of velocity W{t, x) 
According to \X] W{t, x) solves the PDE: 

fdW , , dW\ d'^W 

1- rvit\T. ^ //. — 



We consider W{x,t) = W{uj{t)x - \{t)) =: W{f). Thus, we have: 

dW 

— = i^{t)x-X{t))W'{f), 

Using the specific form of p{x,t) (Uj), we obtain from P^ : 

(16) {a{t) + g{t)x){{Lu{t) - a{t)cj{t))x - A(i))W^'(/) - ^lu;' {t)W" {f) . 
Equation (|16p implies: 

(17) w"2(t)(cr(t) + .g(i)a;)((w(i) - a{t)uj{t))x - A(i)) = {cj{t)x - X{t)f 

Under this condition we find (jj(t) and X{t). From (fT7|) we obtain: 

(18) 5(i)(c:;(t)-a(t)cc>(i))=c^4(^)^ 

(19) - A(i)5(i) + a{t){Lb(t) - a(i)w(t)) = ~2uj^ {t) X{t) , 

(20) - A(i)a(i) = X^{t)u^{t). 

^1/2 

In this system uj(t) and A(t) are unknown and cr(i) = — . q(t) = 

^ ^ ^ ^ ^ ^ VWT2BtTC^ 

WI^T^bT^. (^^^ ® and d). 
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Using HH) and ^, we find 

(21) Am.-^.(,).-i<!^?!±M±^.(,, 

g{t) po t + B 

Let us denote s{t) — iLj^^{t). Thus, from (fTS]) we get 



s(t) + 3a(i)s(i) = -3g"^(t). 
It can be readily concluded that 

(22) ^2^'(^+ :gL (f(^) + C3)-V3 

where 
(23) 
F(i) ^ -^i^i(t + S) Vi" + 2St + C2 + ^Ki{C2 - B^) \n\t + B+ ^t^ + 2Bt + C2I 

and 

^i = -3^^, C3 = ^(0) - F(0). 
From (1211) we find 

(24) \{t)^K2{F{t)+C:,r^'\ 

. ^ UCl'^K^ 
where A 2 — = . 

Finally wc show that p^ - ((20)) are compatible. Substituting uj[t) and \{t) (see 
(1^ and dMl)) in (QOl), we obtain 

(25) f (f (t) + C,)-'" F'it)^ = iff (F(^) + C,)-'/' ^^ 



(here g(t) = = VW+2Bt + C^). 
Further, ^ gets 



^ ' 2 ^^ ' 2 g(t) 2 ^ ' q{t) q{t) 

Therefore, it is easy to see that (|25p holds identically. 

For u}{t) and A(i) defined from HT]) we have ^J.W"{f) = /^^^'(f). Integrating 
gives 



W{f)^W'{0)Jexpf^A df, 
where f{x, t) = uj{t)x — X{t). 
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2. Case 1^0, m^O 

We will find the velocity vector and the density in the following form: 

(a(t)x + 7(i)v\ 
at)x + (3{t)y , P(i, x) = <j{t) + g{t)x - <j{t) + ^{a{t) + [3{t))x. 
W{x,y,t) ) 

f dW dW 

Here the divergency is ait) + /3(i), the vorticity is — — , — - — , f(t) + ^(t) 

\ ay ax 

Firstly we find the functions p{x,t), a{t), /3{t), j{t) and £,{t). 
The conservation of mass ([T]) yields: 

a + gx + {ax + jy)g + {a + gx){a + /?) = 0. 
Therefore 



(27) g + ag + g{a + P) = 0, 

(28) cr + (7(a + /3) = 0, 

(29) 57 - 0. 

From the vorticity conservation equation ([7]) we have the following equations: 

\dx'^dy dy'^ J 
\otox ox ox'' Oy OxOy J 

p{i-i+[a + PM - 7) - K« + /3) + "^(e - 7)) + 

(32) +^(^^x + (iy+{P + m){^x + (iy) + {^-l) {ax + 7j/)) = 0. 
From (I32I) we get 

(33) g{2i - 7 + (a + /3 + m)(2^ - 7) - 2la - l(i) = 0, 



T 
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(34) g0 + /32 + m/? + (e - /)7) = 0, 

(35) a(C - 7 + (a + /3 + m){S, - -f) - la - 1/3) = 0. 

Below we treat particular cases separately. 

1. g{t) ^ 0, a{t) ^ 0. It follows from ^ that -/(t) = 0. In this case instead of 
system ([55)) - (l55|) we have: 

(36) 2^ + 2(a + /3)^-2/a + 2mC-//3 = 0, 

(37) P + (3^ + mP ^ 0, 

(38) ^ + (a + P)^ ~ la + mC - IP = 0. 

From (|36p and ([55]) we get immediately that P{t) = 0. To find functions a{t), 
^(t) and o'(i) we have the system of differential equations obtained from ((27|) . (|28|) 
and jaSl): 

d + 2a2 = 0, 
tr + acr = 0, 



i+{a + m)^ -la = 0. 



Integrating gives 



(39) a{t) = ^ 



2t + Ki 



where Ki = {a{0))-\ 



(40) a{t) = ^' 



VWTkTi' 



where K2 = <7{0)./\Kr\ - o-(0)v/|KO)M 
(41) e(i) = C{t) 



V\^t + Ki\ 



where C(i) = K3 + I f , dr. 
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If TO = 0, then we can integrate (J4T|) : 

(42) ^(i) = , ^^ + 1, 

where K^ = (^(0) - 1)^/\K:\ = (^(0) - VK^W)"'!- 

Thus, for g(t) ^ 0, a{t) ^ and m = we find the following solution for p{x,t) 
and v(t, x, y): 

,.r,^ / X Cr(0) fln X 

(43) p(i,a;)- ^ ^ ' '" 



^(0)V|2t+(«(0))-i| C/2t+(a(0))-i' 
(44) ^^="W'^ + ^(^)^=2^T(I(0JFT' 

(45) V. . e(t). + my - (i + M^l=ilM£))=:f!^ .. 

V v/|2i+("(o))-M y 

2. .g(i) ^ 0, cr(t) = 0. From ^ we get -/{t) = 0. From ^ we find /3(t) as 
follows: 



'0, /3(0) = 0; 

(46) /3(0 = <( ;^?!?=^' /3(0)^0to^0; 



6(0) 



where Ci : , , n ■ 

m + 6(0) 

It is easy to see from (|27|) - (|291) and (l33l) - (|35)) that functions a{t) and ^{t) solve 
the following system: 

(47) i+ia + 0)^-la + m^-'^f3 = 0, 



(48) a + l3 + a{a + l3) + {a + l3f = 0. 

Further, assume to = 0. Using ([M|) and (|T7)) - (P51) . we obtain 



(49) ^ + ^a + P)^~la-^-p = 0, 



(50) a + 2a^ + 3af3 = 0. 
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If a(0) =0 then the solution of ((50|) is zero identically. The function £^{t) can be 
found from (l49l). Therefore 



(51) m= ''^''' 



2(i + C2) ' 
where 

(52) C2 = (/3(0))-\ C3 = 2(/3(0))-ie(0). 

Thus, in this case we obtain 



(53) p{t,x)- ^°^ 



U{t + C2) ' 

('') "^^'' "^^^(TTc^^+rrc^^' 

where the constants C2 C3 are determined early (see ([5^ ). 
If a(0) 7^ 0, we denote A(t) = (a(i))~^ and we get from ([501 

i(t) - 3/3(i)A(i) -2 = 0. 
Therefore, we can find 

where 



(56) '' 



c.^m)r\ c. = M5)^±^. 



ci 



Then we can find from 



where 



(58) C, = {m-\)\j\CACl-l\. 



l_ 
2' 
Thus, if g{t) ^ 0, cr(i) = 0, m = and a(0) 7^ we have 



.r„N ., X PO Ci{t + C2) 

(59) p{t,x)- 



u {t + C2){Ci{t + C2Y-iy 
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(60) Vi 



it + C2)iCiit + 02^-1)-' 



where constants C2, C4 C5 are determined in ([5S)) and (|55|) . 

2.1. The case W{x,y,t) — uji{t)x + uj2{t)y + i^sit). Early we find the first and the 
second component of velocity vector. In this section we consider the special form 
of the third component: 

W{x,y,t) =LUi{t)x + uj2{t)y + U3{t). 
Here ^ and ^ imply: 

(62) {a + gx){Lb2+-/uji+ I3uj2) = 0, 

(63) 

g {ujix + LU2y + tOs + {ax + "fy)uji + {£,x + I3y)uj2) + (ct + gx){LUi + auji + ^^2) = 

We treat particular cases separately: 

1. g{t) ^ 0, a{t) ^ 0. From §7^ and (USD we obtain: 

(64) LUi + auJi + £,UJ2 — 0, 

(65) UJ2 + l3uj2 + ju;i = 0, 

(66) LU3 = 0, 

Early we proved that in the case g{t) ^ 0, a{t) ^ we have /3(i) = and j{t) = 
(see dm) and gSl)). Thus, from ^ and fMl) one gets 

Lu,{t)=uj,{0), i = 2,3. 
Further, we can find from (|64p 



iOi{t)=K{t)- ^ 



jr(^) ^ M^e-"* / ^^f\_dt-l^WTK;\+K,. 
m J ./2t + K, 
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where Ki ^ {a{Oj)-\ 
If TO = 0, then 



K{t) = -iJ2iO)K^t - ^lLU2m2t + Ki I'/' + Ke 



where K^ = (C(0) - 1)^^71, K^ - -^(0) + -?L^2(0)i^i)v^l^- 

2. g(t) ^ 0, a{t) = 0. Then functions uJi{t), i — 1,2,3 solves the system ([M)) - 
(IM)) . It is obvious that tj3(i) = ^3(0). Functions u!i{t) and u!2{t) wc find for case 

TO = 0. We have from dM]) and (gl) 7(i) = 0, /3(i) = t-tttttttt- Thus, it is easy 

t+(/3(0)) 

to see that 

where C2 = (/3(0))^i. 

Further, we can find uji{t). Namely: 



It + Cs 



a). If q(0) = 0, then a{t) = and ^(i) = % , where C3 = 2C2^(0), 

2(i + (72) 

C2 = (/3(0))~i (see dSUand ([52])). Therefore, we obtain from §^: 

(67) c^i(t) = -(i/cc;2(0)i + Ce In |i + C2I + C7), 

where Ce = MO^m " ^)/3"H0), C^t = -C6ln|C2| - c^i(O). 
b). If a(0) 7^ 0, then it follows from ^ that 

|t+C2| 



(68) - / a(r) dr = In ■ 



V\c7it + W' 

u 

where C2 = (j3{0))-\ C4 



(a(0))-i+C2 



Using dHI), jSH) and ([Ml), we obtain: 



V|C4(i + C2)2 

where 



Kit) = Cs — -64 arcsm 



2 .J\Cl\{t + C2) 2(t + C2)2 

Cs - y'|C4C|-l|/3(0)(c^i(0) + /3(0)) + i^arcsin ■ '^^^^ 



2 ^Cy 

We find all components of the velocity field ((26)) . where M^(j:,j/,i) = uoi{t)x+ 
+L02{i)y + ^z{t)- Therefore, we are ready to study the problem of the gradient 
catastrophe: 
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1. If a(0) + 0, /3(0) = and 7(0) = then 

o-(O) po 2^ 



p(t,a;) 






Vl = „, , ^^ , V2=N + 



^3 - f ^^^^=P^ - ^^-^(0)(2^ + K,^ . + .,(0), + .3(0); 

V VI2i + ifil 3 y 

It is obvious that for such velocity field the gradient catastrophe takes place in 

the time t = -— if q:(0) < 0. 

2a(0) ^ ^ 

2. If cr(0) == 0, a{0) ^ 0, /3(0) ^ and 7(0) = then 

., . pox 

p[t,x) = 



U{t + C2) ^ 



"^ = '' "^ = ^(^TC2)" + ^TT^^' 

V3 = -(^;cc>2(0)t + Ce In |t + C2I + C7)x + ;J:^y + a>3(0); 

Z t + C/2 

In this case the necessary condition of gradient catastrophe is /3(0) < and its 
time is equal t = — (/3(0))~^. 

3. If ct(0) = 0, a(0) ^ 0, /3(0) 7^ and 7(0) = then 

,, , PO C4(t + C2)2 

P(i, a;) - 



Vl 



C/(t + C2)(C4(i + C2)2-l)" 



(t + C2)(C4(t + C2)2-l) 



""' V|C4(t + C2)2-l| r t + C2^' 



V3 = (Cs(t + C2) G4t + ^2) 



+ ^y + C.3(0). 
I + 02 

To solve the problem of the gradient catastrophe we find the initial values of 
velocity such as the function f(t) = {C4(t + 6*2)^ — l)(i + C2) has positive zeroes. 

Zeroes of this functions are t = Ti = — C2 and t = T2^3 = ± — C2. T2 > and 

VC4 
r3 > if: 
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■« + i>o 

a(0) + ^ ^ "' 

(/3(0))-3/2 1 



v/(a(o))-i + (/3(o))-i my 

Thus, we obtain: 

1). If /9(0) > and q(0) > — /3(0) then gradient catastrophe doesn't appear; 

2). If /?(0) > and a(0) < — /3(0) then gradient catastrophe appears in the time 



3). If /3(0) < and < a(0) < — /3(0) then gradient catastrophe appears in the 
timeT = -(/3(0))-i; 

4). If /3(0) < and a(0) > — /3(0) then gradient of velocity turns in infinity in 
two points: 



t^T.^ 1 i / «(0) 



/?(0) y^ a(0)+/3(0) ) 

m < T2); 

5). If /3(0) < a(0) < then gradient of velocity turns in infinity in three 
points: 



t = T =^( / "(0) 1 

' /?(0) lV«(0) + /3(0) , 



t = T2 = -{m) 



-1 



,^y^_ 1 ( I aiO) 



m I v«(o)+/9(o) 



(Ti < T2 < T3). 



2.2. The case W{t,x,y) = wi(t)x2 +tj2(i)a;2; + W3(02/^ + ^i(0^+ +A2(i)y + A3(i). 
In this section the third component of velocity vector is function W{t, x, y) ~ 
^LJi{t)x^+Lj2it)xy + Lj3{t)y^ + \i{t)x + \2{t)y + \3{t). We get from pl| and (PT|) : 

(69) 5(^1 + 2auji + C^2) = 0, 

(70) g{uj2 + ia + (3)uj2 + 2"fu;i + 2^u;3) =0, 



14 ANASTASYA KORSHUNOVA 

(71) (ct + g) (w3 + 2/3w3 + 7CJ2) = 0, 

(72) (25 + a)(Ai+aAi+eA2)=0, 

(73) cr(w2 +a + (iuj2 + 27W1 + 2^W3) + (.g + cr)(A2 + /SAa +7^1) = 0, 

(74) gA3 + cr(Ai +q;Ai+^A2) =0. 

To solve this system consider several cases: 

1. g{t) ^ 0, a{t) ^ 0. Then we get from dH]) and (gS]) /3(i) = and 7(t) = 0. 
Therefore, we solve the following system: 

(75) wi + 2awi + iuj2 = 0, 

(76) Ai+aAi+^Aa = 0, 

(77) 0^2 + au}2 + 2^W3 = 0, 

(78) ^3 = 0, 

(79) \ = Q, z = 2,3. 

It easily follows from ^-^ that cj3(t) ee ^3(0), A,(i) ee A,(0), i = 2,3. 
Further, we find from ([76|) and ((77|) : 

We have from dSl) 



(80) / a{T) dr = In v/|2f + -ft:i| 

Using dH]) and dM]), we can solve ([75]) and (|75)) : 



X^{t) = -X2{0) Ui2t + K,)+ ^'' + ^' 



3' ^\2t + K,\l' 
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-t^^\Zt + Kl' 



2{2t + Ki) 3 "^ ' ' 2t + ii:i' 

where 

Ai(0) , ZXi 



if7--^i^ 



A2(0) 3 



2. 5(t) ^ 0, a{t) = 0. We get from ^ and d^l) 7(i) = and /3(t) 
— j-. Thus, we obtain the system: 



(81) uji + 2auji + ^LU2 = 0, 

(82) Ai+aAi+eA2 = 0, 

(83) c02 + (a + /3)tJ2 + 2^t^3 = 0, 

(84) A2 + (3X2 = 0, 

(85) W3 + 2/3w3 = 0, 

(86) A3 = 0. 

Firstly, we get from ^ that A3(t) = A3(0). 

We have from (l46ll 



t 

PiT)dT = ln\t + C2\ 


Therefore, it can be easily calculated that 



(87) A2(t) = A2(0) ^ 



W3(i)=W3(0) 



t+C2 

1 

{t + C2)2 
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Further, we have two cases: 

a). If a(0) = then from ((50l) a{t) = 0. For this reason we have the foUowmg 
system: 

(89) Ai+eA2 = 0, 

(90) wi + Cuj2 = 0, 

(91) LJ2 + PUJ2 + 2^UJ3 = 0, 

where m = ^^^f^ (^^^^ <EID)- 
Let us denote 



where C2 = im)'\ C3 = 2(/3(0))-i^(0) (see dSlD). 
Using dig), find the solution of §^ and dHT]): 

(93) Ai(t) = -M^(/(i)-C9), 



(94) uj2{t) = -uj3{0)(I{t)-Cio), 

where constants Cg and Cio are determined from the initial data, namely: 

Finally, we solve ((90)) : 



(95) ^i(i) = ^f]?(i), 

where 



nl(u) = -C2l{l - I) In^ \u\ + Pu\n \u\ + Cn In |u| + {ICio - I'^h + — + C13, 
2 u 

(here m = < + C2). 
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Cn=C2Cioa-l) + l{lC2~C3), 



C\2^C2il-l)ilC2~C3), 



Ci3 - ^^ - IC21{1 ~ Oln' IC2I - (/'C2 +C„)ln|C2| - {IC^o ~ P)C2 - ^. 

L03(U) Z G2 

b). a(0) ^ 0, then a(i) == {{t + C2){Ci{t + €2)^ - l))-^ (see ([55])) and dM]). We 
have: 



/ a{T) dr — hi ■ 



t + C2 



Ai(i), wi(i) and u}2{t) solve the system (|8T1) - ((83)) . 

The equation (fSTI) is the same as ([5^ . Therefore, we find: 



Xiit)=L{t) 



t + C2 



V|C4(i + C2)2-l| 



e', 



where 



1 



1 



L(i) = C14 Ca arcsin — 

2 " y/\C4\{t + C2) 



Ci4 = VI^^^2-l| (2C2Ai(0) - 1) + Iq arcsin ^^. 

^-2 ^ Vl^4|C2 



We find from 



v/|C4(i + C2)2-l| 

where the function ^2{t) can be found from: 



(96) 



n'^{t) = -{iV\C4{t + C2)^-i\ + 2C5) 



[t + C2f 



Let us denote u ^ t + C2, p{u) — y/C^u'^ ~ 1. Further, we integrate 
we obtain: 



Thus, 



n2{u) - C^3(0) I — + M^1!^L^ _ l^\C^\\ll I V^U+v/|C4li2-l| I+C15 

where 



Ws^Uj V O2 
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Finally, we have from ([82]) 



wi{t)^Wi{t)- 



where Wi [t) can be found from 



w[{t) - -{W\Ci{t + C2Y - 1| + 2C5)- 



^2{t) 



{t + C2Y ■ 

It can be easy calculated that 

Wiiu) = C4 ( piu) + M + arctan^-— + In lul ) + ^/cl\n \\fclu +_p(w)| + 

1 — p(u) 1 ^, , 

where 

J(u) = / ^ -^^ \^ -ln\WC4U + p(u)\du. 

J u^ 

Let us remark that conditions of the gradient catastrophe are the same as for 
linear function W{t, x, y). 

2.3. The case W{t,x,y) = ui{t)x + ^3(^)2/ + W3(i)z + Ui{t). We wiU find the 
velocity vector and the density 



t{t)x + ^{t)y\ 



Po, 



(97) v^\mx + my\. p{t,x)=ait)+git)x^a{t) + i^{a{t)+P{t))x 



W{t,x,y,z) J 
Here the divergency is 



U 






dz 
the vorticity is 

( dW dW ' ^ 



As early, we write the conservation of mass and the vorticity conservation equa- 
tion for our class of solution. Equation ([T]) gives 

dW 

a + gx + {ax + jy)g + {a + gx){a + /3 + -jj—) = 0. 

Therefore 



dW 

(98) g + ag + g{a + (3+—-) = 0, 

oz 
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dW 
(99) & + a{a + l3+—) = 0, 



(100) 57 = 0. 
From ([7|) we have the foUowing equations: 

fd^W dW , ^d'^W BW ,^ ^ ,d'^W\ dW dW 

\otoy ox oxoy oy oy^ J oy oz 

92^ / d^w d^W d^W 

(101) +W--— = ^l' 



dydz \dx^dy dy^ dyd^zj 

\otax ox ox'^ oy oxoy oxoz ' 

(102) = -Ai 



\ Qx"^ dxdy"^ dxdz 

P (e - 7 + (« + /3)(e - 7) - Ua + /3) + m(^ - 7)) + 

(103) +^[Cx + Py+iP + m) i^x + Py) + (C - [ax + 7^)) = 0. 

We consider the case W{t,x,y) = uji[t)x + uj2{t)y + cJ3{t)z + LJ4{t). For such 
function (fTUT|) and (fTU^ imply: 

(104) (cr + gx){uj2 + 7CJ1 + /?a;2 + (^2^3) = 0, 

5(cJia; + uj2y + UJ3 + (ax + 7y)wi + {^x + Py)uj2 + ujz{u}i{t)x + L02{t)y+ 

(105) + ujy,{t)z + LOi{t))) + (cr + gx){u}i + auoi + ^0^2 + tJi^a) = 

(|103|) is the same as ((32)) . Therefore, from (|103p we obtain the system (|33|) - (|35]) . 
Further, we find the functions a{t), a{t), P(t), 7(i), S.{t), uJi{t), ti;2(i), i-^sit) and 

W4(i). 

We can integrate the system (|104|) - (|105p in the case g{t) ^ 0, (7{t) ^ 0. We have 
from PHOI 7 EE 0. Then we may conclude from ^ and ^ that /3 ee 0. 

Thus, we have the system, which follows from (|104p - (|105p and facts 7 = 0, /? = 0: 
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(106) uji + {a + cJ3)uJi + C^2 = Oi 

(107) ^2+^3^2=0, 



(108) W3+u;| = 0, 



(109) LU4 + UJ3LU4 = 0, 

We solve the system ()107p - (|109p and obtain: 



(110) iU3{t) = ^ 



t + cs' 
where C3 = (a;3(0))~^, and 



(111) UJ^{t) 



t + C3' 

where q = Wi(0)c3, j = 2,4. Here 



(see (IW1) V Therefore, we get from 

(112) a + 2a^ + Lu^a = 0. 

(113) <T + cr(a + W3) = 0, 

Let us denote A{t) — {a{t))^^ and obtain from (|112[) 

A - UJ3A -2 = 0. 

Now it can be calculated that: 



(114) ait)^iA{t))-' * + ''' 



t2 + 2c3t + C5 ' 

where C5 — (q;(0))^^C3. 
We have from pi4p : 

t 
(115) a{T)dT = -lii\t^ + 2c3t + C5\. 
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Using pTS]) . we solve (|TT3l) : 

(116) (7(t) = , ^^ 

where cq = cr(0)c3y^. 

Further, we get from ((35|) 



C + a^ - ^a = 0. 
Consequently, 



(117) CW=/+ "' 



Vt^ + 2c3t + C5 ' 

where cy = (^(0) - 0\/c5- 

Finally, we find cji(i) from (|106p 



(118) a;i(i) = ^^*^ 



(t + C3)Vi^ + 2c3t + C5 ' 

where 



X(t) = -c2C7t - -C2I ((t + C3)g(i) + (C5 - cl) In |f + C3 + g(t)|) + ci, 



(here g(t) = Vt^ + 2c3i + C5), 

Cl = Wi(0)c3^/ci'+ -C2Z(c3yci'+ (C5 - C3 

Thus, in this case we obtain: 



(119) pit x) = ^6 P^{t + C3) 

^ ^ ^^ ' ^ (t + C3)Vt^ + 2c3t + C5 + C/(i2 + 2c3t + C5) '^' 

(120) VI = , ^;^''' X, 



(121) V2=(/ + ^= ""^ )x. 

Kit) Cl 1 C4 

(122) V3 = - , ,,.y ^ , =x- + -^y+— z-' 



(i + C3)V<2 + 2c3t + C5 t + C3 t + C3 t + C3' 

where the function -?C(i) and the constants q, i = 2, ..., 7 are determined early. 

To study whether the gradient catastrophe arises, we find the positive zeroes of 
/(t) = = (t + C3)(i^ + 2c3i + C5). Let us denote B — C3, —BC — C5. Thus we obtain 
the problem just considered in. Section 1. 
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